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Abstract 

This paper addresses the issue of homogenization of linear divergence 
form parabolic operators in situations where no ergodicity and no scale sep¬ 
aration in time or space are available. Namely, we consider divergence form 
linear parabolic operators in C M n with x (0, T))-coefficients. It 

appears that the inverse operator maps the unit ball of L 2 (fl x (0,T)) into 
a space of functions which at small (time and space) scales are close in H 1 - 
norrn to a functional space of dimension n. It follows that once one has 
solved these equations at least n-times it is possible to homogenize them 
both in space and in time, reducing the number of operations counts neces¬ 
sary to obtain further solutions. In practice we show that under a Cordes 
type condition that the first order time derivatives and second order space 
derivatives of the solution of these operators with respect to harmonic co¬ 
ordinates are in L 2 (instead of H -1 with Euclidean coordinates). If the 
medium is time independent then it is sufficient to solve n times the associ¬ 
ated elliptic equation in order to homogenize the parabolic equation. 
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1 Introduction and main results 


Let 17 be a bounded and convex domain of class C 2 of W 1 . Let T > 0. 
Consider the following parabolic PDE 

f dtu = div (a(x, t)Vu(:c, t)) + g in f2x(0, T) 

| u(x, t) = 0 for (x,t) G (<90 x (0, T)) U (O U {t = 0}). 

Write Ot := O x (0, T). g is a function in L 2 (£1t)- (x,t) —> a(x,t) is a 
mapping from VLt into the space of symmetric positive definite matrices 
with entries in L°°(Qt)- Assume a to be uniformly elliptic on the closure 
of ftp. This paper addresses the issue of the homogenization of CD in 
space and time in situations where scale separation and ergodicity at small 
scales are not available (see [13, ED and (3] for an introduction to classical 
homogenization theory). For that purpose, we will introduce in subsec¬ 
tion rm theorems establishing under Cordes type conditions the increase 
of regularity of solutions of CD when derivatives are taken with respect 
to harmonic coordinates instead of Euclidean coordinates. In subsections 
oo these results will be used to homogenize CD in space and in time. 
More precisely, assume a to be written on a fine tessellation with N degrees 
of freedom. If a is time independent, then by solving n-times an elliptic 
boundary value-problem associated to (11.11) (at a cost of 0(IV(ln N) n+3 ) 
operations using the Hierarchical matrix method ED) it is possible approx¬ 
imate the solutions of CD by solving an homogenized operator with N a 
degrees of freedom (a < 1, a = 0.2 for instance) or with a fixed number 
M of degrees of freedom (numerical experiments given at the end of this 
paper have been conducted with N = 16641 and M = 9), this problem is of 
practical importance for oil extraction and reservoir modeling in geophysics. 
If a is also characterized by a continuum of time scales, then the method 
presented here does not reduce the number of operation counts necessary to 
solve (11.11) only one time. However if one needs to solve (11.11) K (K > n) 
times (with different right hand sides) then by solving (11.11 ) n times it is 
possible to obtain an approximation of the solutions of CD by solving an 
homogenized (in space and time) parabolic equation written on a coarse 
tessellation with coarse time steps. 

1.1 Compensation phenomenon 

Let F be the solution of the following parabolic equation 

dtF = div (a(x, t)VF(x, f)) in 0^ 

< F(x,t) = x for (x,t) G (dfl x (0, T)) (1.2) 

k div (a(x, 0)VF(x, 0)) = 0 in fl. 
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ByO) we mean that F := (F \...., F n ) is a n-dimensional vector field such 
that each of its entries satisfies 


dfFi = div ( a(x,t)VFi(x,t )) in f It 

< Fi(x,t) = Xi for ( x,t ) G (dft x (0, T)) (1-3) 

div (a(x, 0)VFi(x, 0)) =0 in ft. 


Observe that if a is time independent then F is the solution of an elliptic 
boundary value problem. 

Definition 1.1. Write 

a := 4 V FaVF . (1.4) 

Write /3 a the Cordes parameter associated to a defined by 

/ (Trace[cr]) 2 \ 

H, := esssup (l , t)Eni . [n - )- (1.5) 

Observe that since 


At = esssup (a . it)eaT 


n — 


(E"=l \,cr(x,t))‘ 

Er=i a 2 


i,cr(x,t) 


( 1 . 6 ) 


where (A^m) denotes the eigenvalues of M, /3 a is a measure of the anisotropy 
Of (7. 


1.1.1 Time independent medium. 

In this subsection we assume that a does not depend on time t. Write 
for p > 2, Wyf ( D for Dirichlet boundary condition) the Banach space 
W“jf{ft) O l / l / Q iP (0). Equip W“jf{ft) with the norm 


^1 W 2 a p m 


Equip the space L p (0, T, W^ p (ft)) with the norm 

r T 


/n 


Zl( did J v ) 2 ) 


(1.7) 


h3 


l Lp(o,T,wi; p {n)) 


f f ( y^(djdj v ) 2 ) 2 dx dt. ( 1 . 8 ) 

Jo Jn 


Theorem 1.1. Assume that dta = 0, g £ L 2 (ftx), ft is convex, (3 a < 1 and 
(Traced) 4 £ L°°(ft) then uoF~ l £ L 2 (0, T, W^ 2 (ft)) and 




C 


\l 2 (o,t,w^ 2 (U)) — 


l-(33 


— \\9\\l 2 (Q t )- 


(1.9) 
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Remark 1.1. The constant C can be written 


C 


C n 

(Amin (a)) ^ 


(Trace[ ct]) 4 


Through this paper, we write 


A min (a) := inf inf t l.a(x,t).l. (1.10) 

(x,t)eCi T ieK",|/|=i 


Remark 1.2. According to theorem ll ^ although the second order derivatives 
of u with respect to Euclidean coordinates are only in L 2 ( 0, T, H~ l ( fl)), they 
are in L 2 (Qt) with respect to harmonic coordinates. 

Remark 1.3. Observe that if a is time independent then F and a are time 
independent and F is the solution of the following elliptic problem: 


divaVF = 0 in SI 

< ( 1 . 11 ) 

^F(.t) = x on <90. 

In dimension one F is trivially an homeomorphism. In dimension 2 this 
property follows from topological constraints ;2] (even with a t .j e L°°(0)), 
[B] (one can also deduce from P that for n = 2, if a is smooth then the 
conditions (3 a < 1 and (Trace[<r]) _1 e L°°( O) are satisfied). In dimension 
three and higher F can be non-bijective even if a is smooth, we refer to P 
and m, however in dimension 3 the assumption (Trace[<r])4 1 g L°°(Qt) 
implies that F is an homeomorphism. If n > 4 we need to assume that F is 
an homeomorphism to prove the theorem. 

Remark 1.4. In fact the condition (Trace(cr)) _1 6 L p (£It) for p < oo de¬ 
pending on n is sufficient to obtain theorem 11 .H and the following compen¬ 
sation theorems. For the sake of clarity this paper has been restricted to 
(Trace(cr))^ 1 € L 00 (flr). 

Remark 1.5. Write 

Amax( cr ) 0 \ 

p a : = esssup QT --—. (1.12) 

AminW ) 

It is easy to check that p a is bounded by an increasing function of (1 — (3a)- 1 
and in dimension two (3 a < 1 is equivalent to p a < oo. 

Remark 1.6. Theorem 11.11 has been called compensation phenomenon be¬ 
cause the composition by F increases the regularity of u E T 2 (0, T, Hq (fl)). 
The choice of this name has been motivated by F. Murat and L. Tartar’s 
work on H-convergence ESI which is also based on on a regularization prop¬ 
erty called compensated compactness or div-curl lemma introduced in the 
70’s by Murat and Tartar EH> EH (we also refer to E53 for refinements of 
the div-curl lemma). 
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Figure 1: Site Percolation 


The compensation phenomena presented in this subsection can be ob¬ 
served numerically. In figure U the value of a is set to be equal to 1 or 100 
with probability 1/2 on each triangle of a fine mesh characterized by 16641 
nodes and 32768 triangles. (EH has been solved numerically on that mesh 
with g = 1. u, uoF~ 1 , d x u and dx(uoF -1 ) have been plotted at time t = 1 
in figure El 
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0 . 015 *. 



(c) d x u. (d) dx(uoF 1 ). 

Figure 2: u, u o F~ 1 , d x u and dx{u o F -1 ) at time t = 1 for the time independent 
site percolating medium. 
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In situations where g G L°°(0, T, L 2 (Q)), dtg G L 2 (0, T, 17 _1 (n)) or 5 G 
L p (Q,t) with p > 2, one can obtain a higher regularity for u o I 7-1 . This is 
the object of the following theorems. 

Theorem 1.2. Assume that fl is convex, g G L°°(0, T, L 2 (fl)), 

G L 2 (0, T, i7 - 1 (hl)) ; c^a = 0, (3 a < 1 and (Trace[<r]) 4 -1 g L°°(TIt) then 
for all t G [0,T], no F _1 (., f) G W^ 2 (Q) and 


\uoF (.,t)|| 




c 


i-fl 


(||£'IL°°(0,tt 2 (^)) + H <9 *5 , llL2(0,T,^ 1 (n)) 


(1.13) 


Remark 1.7. The constant C* can be written 


C 


Cn,n 

(Amin(O')) 4 


(Trace[cr])4 1 


it 1 + 


L °°^ y A min (a) 


1 

2 . 


Theorem 1.3. Assume that fl is convex, g(., 0) G L 2 (Q), 
dtg G L 2 (0, T, i7 _1 (n)), g G L p (Qt), dta = 0, (3 a < 1 and (Trace[cr]) 4 _1 g 
L 00 (ri 7 ’) t/ien t/iere exists a real number p$ > 2 depending only on n,Q and 
(3 a such that for each p such that 2 < p < po one has 


LP(o,T,wi p (d)) — i (ll5 , llix(n T ) 

23 l-(33 (1.14) 

+ ||fi , (-lO)|| L 2 ( . f ^ + ||^|| i 2( 0)TH -l( n ))). 


Remark 1.8. The constant C can written 


C = 


a 


n,Q,p 


('^min(flO) 4 


w (Trace[cr])4 


4-H 


)( 1 + 


L °°^ K A min (a) 


)*■ 


Write 


Tllcr(a) := sup 

x,yEQ,Xy£y 


\v(x) - v(y )I 


\x - y |" 


(1.15) 


Theorem 1.4. Assume that n < 2, Tl is convex, g(., 0) G L 2 (Q), dtg G 
L 2 (0,T, 17 _1 (n)) ; g G LP(fi r ), d t a = 0, (3 a < l, (Trace^])" 1 G L°°(h! T ) 
and g G L 2 [0,T; L p * (Q)\ with 2 < p*. Then there exists p G (2 ,p*] and 
7 (p) > 0 such that 

( f ||V(noF" 1 )(.,t)||^, n) dt ) 2 < ° 1 (Hff||Lp(n r ) 

Jo 1 -(33 

+ || 5 (-) 0 )|| L 2 (q) + II^ 5 , IIl 2 (o,t,h- 1 (Q)))- 

(1.16) 
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Remark 1.9. The constant C in (Qj£D depends on n. p, 11 , A min (a) and 
||(Trace(cr)) _1 || L00 ^^. It is easy to check that if n = 1 then the theorem is 
valid with 7 = 1 / 2 . 

In the following theorems Ll is not assumed to be convex. 

Theorem 1.5. Assume n > 2 and dta = 0. Let p > 2. There exist a 
constant C* = C*(n,dLl) > 0 a real number 7 > 0 depending only on n,Ll 
and p such that if (3 a < C* then 

(i l|v(“» Oilmen , dt ) i S c(|| 9 ||„ (nT) + list •’°)IL 2 (n) (117) 

+ ll^5llL 2 (0,T,H- 1 (n)))- 

Remark 1.10. The constant C in (firm depends on n, 7 , C*, A m ; n (a) 

and II(Trace {a)) 2 ? || Loo(Qt) . 

It is easy to check that if a = e(x)S(x,t) where e is a time indepen¬ 
dent symmetric uniformly elliptic matrix with L°°(f2) entries and S is a 
regular uniformly positive function then the results given in this sub-section 
and the homogenization schemes of sub-section 11.21 remain valid with the 
time independent harmonic coordinates associated to e, i.e. solution of 
— div eVF = 0. 


1.1.2 Medium with a continuum of time scales. 


In this subsection the entries of a are merely in We need to 

introduce the following Cordes type condition. 

Condition 1.1. We say that condition 1 7. ?l is satisfied if and only if there 
exists 5 E (0, 00 ) and e > 0 such that 


esssup^ T 


S 2 Trace ^cra] + 1 
(5 Trace [u] + l) 2 


1 

< - . 

n + e 


(1.18) 


Write 


Z<j 


esssup QT n 


Trace/acr] 
(Trace [cr]) 2 ’ 


(1.19) 


Observe that z a is a measure of anisotropy of a, in particular 1 < z a < n 
and z a = 1 if o is isotropic. Write 


y a := || Trace[a]|| i cx J(f2T) || (Trace[cr]) 1 || i0O( n T) - ( L20 ) 

Proposition 1.1. If || Trace[cr]||£<x>(Q T ) < 00 and ||(Trace[<7])~ 1 || L00 ^^ < 
00 then condition \l.l\ is satisfied with 


6 :=n||(TVace[a]) 


and with e := 2n J Ja 2 n provided that z a < 1 + —. 

71 


( 1 . 21 ) 








Remark 1 . 11 . Observe that in dimension one z a = 1, thus for n = 1 condi- 
tion ll.ll is satisfied is Trace[<r] E L°°(Qt) and (Trace[cr ]) _1 G 
We have the following theorems 

Theorem 1.6. Assume that 12 is convex, and condition \l.l\ is satisfied then 
u o F~ 1 G L 2 (0,T, W^ 2 (12)), dfiuoF- 1 ) G L 2 (12 T ) and 

\\ uoF ~ l \\ L H 0 ,T,wZ 2 m + ll^( woi? ” 1 )llL 2 (n T ) < C\\g\\ L 2 (nT) (1.22) 

where C depends on 0, n, 6 and e. 

Remark 1.12. According to theorem 11.61 although the second order space 
derivatives and first order time derivatives of u with respect to Euclidean 
coordinates are only in L 2 (0,T, H^ 1 (Q)), they are in L 2 (Qt) with respect 
to harmonic coordinates. 

Similarly we obtain the following theorems in situations where g G 
L p (Qt) with p > 2 . 

Theorem 1.7. Assume that 12 is convex, and condition o is satisfied then 
there exists a number po > 2 depending on n,12,e such that for p G (2,po)> 
u o F~ l G L p (0,T,W% p (n)), dfiuoF- 1 ) G LP(12 r ) and 

11“ ° F \\LP(0,T,wfi p (Sl)) + II dfiuoF 1 )||lp(h t ) < C'|| 5 lUp(f 2 T ) (1-23) 

where C depends on 12, n, 5 and e. 

Theorem 1.8. Assume that 12 is convex, and condition o is satisfied then 
there exists a number ao > 2 depending on n, 12, e such that for a G (0, afi), 
V(u o F~ l ) G L 2 (0, T, C'“(12)) and 

||V {uo F 1 )(.,t)||L 2 (o,T,C“(a)) < C'll5 , llLp(n T ) (1-24) 

where C depends on 12, 5, n, and e. 


These compensation phenomena can be observed numerically. We con¬ 
sider in dimension n = 2 , 


a(x,y,t ) 


1 1.1 + sin(27rx'/ej) 

6 1.1 + sm(27ry'/ei) 


+ sin(4x' 2 y /2 ) + 1) 


(1.25) 


with x' = x + \/2t, y' = y — y/2t, e± = g, e 2 = jg, e 3 = gy, £4 = gg and 
65 = gg. This medium has been plotted in figure 0 at time 0 (observe that 
^max (a)/A m ; n (a) ~ 100 ). 

EU) has been solved numerically on that mesh with g = 1 on the fine 
mesh characterized by 16641 nodes and 32768 triangles. Figure^shows d x u 
and d x (u o F~ 1 ) at time 0.3. 

In figure El and |H1 the value of xq is set to xo := (0.75,-0.25) and the 
curves t —> u{xQ,t),uoF~ l {xQ,t), Vu(xo,t), VuoF~ 1 (xo,t) are plotted from 
t = 0 to t = 0.3 
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Figure 3: a at time 0. 




(a) d x u. (b) d x (uoF *). 

Figure 4: d x u and d x {u o F _1 ) at time t = 0.3 for the Multi-scale Trigonometric 
time dependent Medium. 
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the point is (0.75,-0.25) 



(a) u and u o F 1 . 

Figure 5: £ —> u(x 0 , t),u o F _1 (x 0 , t) from t = 0 to t = 0.3 with x 0 := (0.75, —0.25) 


the point is (0.75,-0.25) 


the point is (0.75,-0.25) 




(a) d x u and d x (u o F 1 ). (b) d y u and d y (uoF : ). 

Figure 6: t —> Vm(x 0 , t), Vu o F^ 1 (x 0 ,t) from t = 0 to t = 0.3 with x 0 : = 
(0.75,-0.25) 
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1.2 Homogenization in space. 


Let Xh be a finite dimensional subspace of iLg(P) fi tL 1 , 0 °(P ) 1 with the 
following approximation property: there exists a constant G'x such that for 


all / E W 2 /(n) 


Jrf^ ||/ - v\\ H i (n) < C x h\\f\\ w 2 , 2 {ny (1.26) 

It is known and easy to check that the set of piecewise linear functions 
on a triangulation of f l satisfies condition ( 061 ) provided that the length 
of the edges of the triangles are bounded by h (Cx hr (11.261) being given by 
the aspect ratio of the triangles). 

For media characterized by a continuum of time scales we will consider 
twice differentiable elements satisfying the following usual inverse inequali¬ 
ties (see section 1.7 of ESI) : for v E Xh, 

IMIu/; 2 (n) — IMItfi(n)- (1-27) 

and 

IMIn:}(n) < CA-/i _1 ||u|| L a (n ). (1.28) 

In this paper we will use splines to ensure that condition (1071) is satisfied 
(observe that it requires the quasi-uniformity of the (coarse) mesh, i.e. a 
bound on the aspect ratio of the (coarse) triangles). 

For t E (0, T ) let us define 


Vh(t) := {ipoF(x,t) : tp € X h }. (1.29) 


Write L 2 (0, T; 77q(P)) the usual Sobolev space associated to the 


norm 


v 


II 2 



(1.30) 


Write Yt the subspace of L 2 (0, T; Hq (P)) such that for each v E Yt and 
t E [0,T], x —> v(x,t) belongs to V/(t). 

Write Uh the solution in Yt of the following system of ordinary differential 
equations: 


I (i/j,dtu h ) L 2 (fy + a[ip,u h ] = (V>,9)L2(fi) for all t E (0 ,T) and ip E V h (t) 
|u ft (a;,0) = 0. 

(1.31) 

Write 

a[v,w\ := / t Vv(x, t)a(x, t)Vw(x, t) dx. (1.32) 

Jn 

1 W 1 ’°° is the usual space of uniformly Lipschitz continuous functions. 
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1.2.1 Time independent domain 

We have the following theorem 

Theorem 1.9. Assume thatdta = 0 , Q is convex, (3 a < 1 and (Trace[cr]) _1 E 
then 

l|( u ~ u h){-, T) || ^2^) + || u — U h\\ L 2 (0,T;H^(CI)) ~ ^hWdW L 2 (Cl T ) ■ (1.33) 

Remark 1.13. The constant C depends on Cx, n , 12, A min (a) and 
11( Trace [cr]) —1 ||x / ° 0 (r 2 r )‘ If n > 5 it also depends on || Trace [<t]|| loo ^^ and if 
n = 1 it also depends on A max (a). 

1.2.2 Medium with a continuum of time scales. 

Theorem 1.10. Assume that 12 is convex, and condition ra is satisfied 
then 


||(m 'ihi)(T)|| i2 ^ + | I'M u h\\L 2 (p,T,Hjj(si)) — Ch\\g\\L 2 (n T )- (1-34) 

Remark 1.14. The constant C depends on Cx, n, 12, 5 and e, A min (a) and 

lmax(®) - 

The system of ordinary differential equations OH is still characterized 
by a continuum of time scales in situations where the entries of a merely 
belong to L°°( 12j’). They need to be discretized (homogenized) in time in 
order to be solved numerically. This will be the object of the next subsection. 
Loosely speaking, although ten is associated to a fine tessellation and fine 
time steps, it is possible to approximate its operator on a coarse tessellation 
with coarse time steps. 

1.3 Homogenization in space and time. 

Let M E N*. Let (t n = n-p)o< n <M be a discretization of [0, T]. Let (y?j) be 
a basis of Xh. Write Zp the subspace of Yp such that w E Zt if and only if 
w can be written 

w{x, t) = Ci(t)ipi(F(x, t)). (1.35) 

i 

and the functions t —> efit) are constants on each intervals (t n ,t n+ i\. Write 
V the subspace of Yp such that its elements can be written 

ip(x, t) = di<pi(F(x, t )). (1.36) 

i 

where the parameters di are constants (on [0, T]). For w E Yp, define w n G V 
by 

w n (x,t ) := y ^Cj(t n )<pi(F(x,t)). (1.37) 
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Write v the solution in Zt of the following system of implicit ordinary 
differential equations (such that v(x, 0) = 0): for n £ {0,..., M — 1} and 

{fP(fn+l )■> v n+l (^n+l)) l 2 (Q) = {fl , (fn) i v nitn)) ^ 2(q) 

ptn -\-1 / 

+ / ({dt'<P(t),v n+ i{t)) l2{q) 

J t n 

~ a[i>(t),v n+1 (t )] + (ip(t),g(t)) L2 ^ dt. 

(1.38) 

The following theorem shows the stability of the implicit scheme (PHI) . 
Theorem 1.11. Let v £ Zt be the solution of arm . We have 

IK t )|Il2(o) + IMIx 2 (o,T,xri(n)) < C , ||s , l|.L2(n T )- (1-39) 

Remark 1.15. The constant C depends on n, fl and A m i n (a). 

The following theorem gives an error bound on the accuracy of time 
discretization scheme PHI) when a does not depend on time. 

Theorem 1.12. Let v £ Zt be the solution of arm and Uh be the solution 
of \1.31\) . Assume that dta = 0. We have 

|| ( u h — t ’)(^)|L2(Q) + ll w /i — ' u llL2(o i T,Hl(n)) ^ C|At| 

( ^ x (1.40) 

(11^11x2(0^-1(0)) + || 5 (- 5 0 )|| L 2 (n) J- 

Remark 1.16. The constant C depends on n, 11 and A m i n (a). 

The following theorem gives an error bound on the accuracy of the time 
discretization scheme PHI) when a has no bounded time derivatives. 

Theorem 1.13. Assume that Ll is convex, and condition ra is satisfied. 
Let v £ Zt be the solution of arm and Uh be the solution of hl.dU) . we 
have 

||(tt/j - u)(T)|| L2( . n , ) + II Uh - ^ C , ^“lls'llL 2 (n T ) (1-41) 

where C depends on Ll, n, 6, e, A m i n (a) and A max (a). 

Remark 1.17. Observe that the accuracy of the time discretization scheme 
PHI) requires that |Ai| << h when a has no bounded time derivatives. 

We refer to section |H] for numerical experiments. 
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1.4 Literature and further remarks. 


For early works on homogenization with random mixing coefficients we refer 
to HU, EH, |58|, EU, EH], EE], EH, El- Papanicolaou and Varadhan m 
have considered a two-component Markov process (x(t),y(t)) where y(t) is 
rapidly varying (and is not assumed to be ergodic in dimension one) and 
enters in the coefficients of the stochastic process driving x(t). They have 
studied the convergence properties of x(t) as the fluctuations of y(t) becomes 
more rapid using the martingale approach to diffusion, developed by Stroock 
and Varadhan m, m, Eg, m- 

The numerical homogenization method implemented in this paper is a 
finite element method. The idea of using oscillating tests functions can 
be back tracked to the work of Murat and Tartar on homogenization and 
H-convergence, we refer in particular to m and EL Those papers also con¬ 
tain convergence proofs for the finite element method in an abstract setting 
for a sequence of //-converging elliptic operators (recall that the frame¬ 
work of H-convergence is independent from ergodicity or scale separation 
assumptions and are based on the compactness of any sequence of solutions 
of — diva e Vu e = g with uniformly bounded and elliptic conductivities a e , 
we also refer to the initial work of Spagnolo m for G-convergence). 

The numerical implementation and practical application of oscillating 
test functions in numerical finite element homogenization have been called 
multi-scale finite element methods and have been studied by several authors 
|H], EH], |M], eu, jdiij. [Mj. EH], E]- The work of Hou and Wu EH has 
been a large source of inspiration in numerical applications (particularly for 
reservoir modeling in geophysics, we refer to EH, EH], El and E23 for recent 
developments) since it was leading to a coarse scale operator while keep¬ 
ing the fine scale structures of the solutions. With the method introduced 
by Hou and Wu, the construction of the base functions is decoupled from 
element to element leading to a scheme adapted to parallel computers. A 
proof of the convergence of the method is given in periodic settings when the 
size of the heterogeneities is smaller than the grid size and an oversampling 
technique is proposed to remove the so called cell resonance error m when 
the size of the heterogeneities is comparable to the grid size. 

Allaire and Brizzi |3| have observed that multiscale finite element method 
with splines would have a higher accuracy and have introduced the compo¬ 
sition rule (we also refer to |H] ) • In j5H], it has been observed that if u is the 
solution of the divergence form elliptic equation 


— div (a(x)Vn(x)) =g in Q 
u = 0 in <9fi. 


(1.42) 
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and F are harmonic coordinates defined by 


div aV-F = 0 in Q 
F(x) = x on dfl. 


(1.43) 


then under the Cordes type condition /3 a < 1 on a given by m, one has 
for some p > 2. 

\\u o F^Wwz.Ptn) < C\\g\\ LP ^y (1-44) 

It has been deduced from this compensation phenomenon that numerical 
homogenization methods based on oscillating finite elements can converge in 
the presence of a continuum of scales if one uses global harmonic coordinates 
to obtain the test functions instead of solutions of a local cell problem m- 
In dimension three and higher it has been known since the work of Fenchenko 
and Khruslov ESI, EH that the homogenization of divergence form elliptic 
operators — diva e Vri e = g can lead to a non local homogenized operator if 
the sequence of matrices a e is uniformly elliptic but with entries uniformly 
bounded only in L 1 (fl). From a numerical point of view this non-local effects 
imply that a nonlocal numerical homogenization method cannot be avoided 
to obtain accuracy. Hence in m , it is shown that the accuracy of local 
methods depend on the aspect ratio of the triangles of the tessellation with 
respect to harmonic coordinates (which is not the case if one uses non local 
finite elements, we refer to m for further discussions on the apparition of 
non local effects in numerical homogenization). Recently Briane has shown 
(20| that this non-local effect is absent in dimension two in the H-convergence 
setting. 

The phenomenon is similar here, however observe that if one has solved 
the initial parabolic equation at least n times and those solutions are (lo¬ 
cally) linearly independent it is also possible to use them as new coordinates 
for numerical homogenization. Observe that in dimension higher than three 
the harmonic coordinates are not always invertible, an idea to bypass this 
difficulty could be either to choose the change of coordinates locally and 
adaptively or to enrich the coordinates by writing down the initial equations 
as degenerate equations in a space of higher dimension [£2], these points have 
not been explored. For divergence form elliptic equations, recall that fast 
methods based on hierarchical matrices 2 are available EH El El EDI EH for 
solving (EH) and eh in 0(iV(ln 1V) TI+3 ) operations (N being the number 
of interior nodes of the fine mesh). 

The issue of numerical homogenization partial differential equations with 
heterogeneous coefficients has received a great deal of attention and many 
methods have been proposed. A few of them are cited below. 

2 As for the fast multipole method and the hierarchical multipole method designed by L. 
Greengard and V. Rokhlin m, these methods are based on the singular value decomposition of 
operators Green’s function. 
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• Multi-scale finite volume methods m- 

• Heterogeneous Multi-scale Methods EH, ESI- 

• Wavelet based homogenization m, m, m, ixg, m, eg. 

• Residual free bubbles methods m- 

• Discontinuous enrichment methods 123 , EH- 

• Partition of Unity Methods m- 

• Energy Minimizing Multi-grid Methods CHI- 

Following the methods of E3, it is possible to implement a hnite-volume 
method based on the compensation theorems given in this paper. The ele¬ 
ments given in this paper contain the fine scale structure of F. as it has been 
done in ESI, it is possible to approximate the initial parabolic operator by 
a homogenized parabolic operator associated to the coarse mesh (the test 
functions in this case would be piecewise linear on the coarse mesh and the 
approximation error associated to the homogenized operator would depend 
on the aspect ratio of the triangles of the coarse mesh in the metric induced 


by F). 


Finally, in this paper a has been assumed to be bounded and uniformly 
elliptic. Without these assumptions the diffusion associated to homogenized 
operator can be anomalously slow m, Eg or fast (super-diffusive) EH • If a 
has an unbounded skew symmetric component, the homogenization of <HU> 
can give rise to a degenerate operator EH- 


2 Proofs 


2.1 Compensation. 

2.1.1 Time independent medium. 

We will need the following lemmas. Let At be the bilinear form on L 2 (0, T ; Hq (O)) 
defined by 



where 


“M(t) : = t Vv(x, t)a(x, t)Vu(x, t) dx. (2.2) 


n 


We write MtM := At[u,u]. 

Lemma 2.1. We have 


IK.r)fe (n) + A T [ul < (2-3) 


(2.3) 
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Proof. Multiplying (iniby u and integrating with respect to time we obtain 
that 

-\\u(., T )\\ 2 L2{ a ) + At[u\ = (u,g) L 2 ^ T y (2.4) 

Using Poincare and Minkowksi inequalities leads us to Q. □ 

Lemma 2.2. Assume dta = 0. We have 

\\ d t u \\ 2 L2(n T ) + a H-’ T )] ^bf L 2 ( n r) - (2- 5 ) 

Proof. Multiplying «ni> by dtu and integrating by parts we obtain that 

\\ d tu{.,t)\\ 2 L2(n) + a[d t u,u} = (d t u,g) L 2 (ny ( 2 . 6 ) 

Observing that 

a[d t u,u} = ^9 t (a[u]) - \ ud t aS7u. (2.7) 

2 2 J n 

we conclude by integration with respect to time and using Minkowski in¬ 
equality. □ 

Lemma 2.3. Assume dta = 0. We have 

\\dtu(., T )\\ 2 L 2 ^ +A T \d t u\ 11^11^(0,r,H-i(n)) + 0 )|lx,2(f2)- 

( 2 . 8 ) 

Proof. We obtain from GH> that 

d^u = div (a(x, t)Vdtu(x, t)) + div (d t a(x, t)S/u(x, t)) + dtg. (2.9) 
Multiplying m by dtu and integrating with respect to time we obtain that 


2 IIj|| L 2 (n) 


+ 


||cW-,0)||^( n ). 


(dta) [ dtu , u] cLt 


( 2 . 10 ) 


We conclude by the H 1 -duality inequality and Minkowski inequality. □ 

We now need a variation of Campanato’s result m on non-divergence 
form elliptic operators. Let us write for a symmetric matrix M, 


Trace (M) 
Trace (* AIM) 


( 2 . 11 ) 
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We consider the following Dirichlet problem: 


L m v = f 


( 2 . 12 ) 


with Lm '■= Yl'i j=\ Mij( x )didj. The following theorems 12.11 and 12.21 are 
straightforward adaptations of theorem 1.2.1 of EH- They are proven in 
m under the assumption that M is bounded and elliptic. It is easy to 
check that the conditions (3m < 1 and vm E L°°(H) are sufficient for the 
validity of those theorems. We refer to m for that adaptation. 

Theorem 2.1. Assume that /3 m < 1, v m £ L°°(Q) and 11 is convex. Then 
if f E L 2 (Q) the Dirichlet problem 12.121) has a unique solution satisfying 


\w 


2,2 (Cl) — 


C 


1-/3 


— \WMf\\L 2 (Q)- 


(2.13) 


M 


Remark 2.1. (3m is the Cordes parameter associated to M. 


Theorem 2.2. Assume that /3 m < 1, vm £ T°°(12) and 11 is convex. Then, 
there exists a real number po > 2 depending only on n, H and (3 m such that 
for each f e L P (H), 2 < p < po the Dirichlet problem \2.12 1) has a unique 
solution satisfying 


W- 


2 ’ p (fl) — 


^AlWMfWmay 

1-/3' 


M 


Let us now prove the compensation theorems. Choose 

a 


M : = 


det(VF)| 2 


o F 


-l 


(2.14) 


(2.15) 


It is easy to check that (3 a < 1 implies that F is an homeomorphism from 
H onto H, thus (12T51) is well defined. Moreover observe that (3m = (3 a and 


\ v M\\l°°(Q. t ) ^ 


C„ 


(A m i n (a)) 2 


g-1| (Trace[< t]) 4 1 ||^ 0 


'(n T )- 


(2.16) 


Fix t e [0, T], Choose 


/: = 


(dtu - g ) 

det(VF)|3 


o F 


-l 


(2.17) 


Observe that by the change of variable y = F(x) one obtains that if dta = 0 
(which implies that F is time independent), dtu E L 2 ( 11) and g(.,t ) E L 2 (H) 
that / E L 2 ( 11) and 


ll/llL 2 (n) < l|c>tu|| L 2(n) + ||5 , llL 2 (a)- ( 2 - 18 ) 
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__ 2 2 

It follows from theorem 12.1 1 that, there exists a unique v £ W r j (Q) satisfying 




< 


C\\i'm\\loo(t 1t ' ) 


(1-P3) 2 

and the following equation 


(11*^11 L 2 (fi) + II^IIl 2 ^))- (2.19) 


dtu(y,t) = ^2 (a(F 1 {y,t),t)) i j d i d j v(y,t)+g(y,t). ( 22 q) 


h3 


We use the notation g := g o F and u := u o F~ l . Using the change of 
variable y = F(x) and using the property div aS7F = 0 when dta e 0 we 
obtain that E3) can be written 


dtu = div (aV(v o F)) + g. (2-21) 

If dtu £ L 2 ( fl) and g(.,t) £ L 2 (Q) we can use the uniqueness property of 
the solution of the divergence form elliptic Dirichlet problem 


div (aVrc) = dtu — g. 


( 2 . 22 ) 


to obtain that v o F = u. Thus using lemma EQ1 we have proven theorem 
o Moreover assume that g £ L 2 (Qt) and dtu £ L 2 { Qt)- It follows that 
for t £ [0, T] — B, g(.,t ) £ L 2 (Q) and dtu(.,t) £ L 2 (Q) where B is a subset 
of [0, T] of O-Lebesgue measure. It follows from the previous arguments that 
on [0,T] — B, u o F~ 1 (.,t) £ W 2 j 2 (iP) and satisfies 


\uoF 1 (.,t)|| 


- 


C\\v M \ 


l°°{ n T ) 


(1-/3I) 2 


(\\d t u(.,t)\\ 2 L 2 ( n ) + 


(•> III 2 (fi))- 

(2.23) 


Integrating with respect to time we obtain that 

u o F~ l £ L 2 ( 0, T, W% 2 (n)) and 


\\ U °F ^L 2 (0,T,W^ 2 (n)) 


< 




(i -m 


■(ll^ u lll 2 (o T ) + lblli 2 (f2 T ))- 


Thus using lemma. 12721 we have obtained theorem II .1 1 

Let us now prove theorem 11.31 Assume that there exists qo > 2 such 
that for 2 < p < qo, dtu £ LP(Qt) and g £ L p {VIt)- Let us now apply 
theorem 12.21 with p < min(po, qo), M given by (12.151) and / given by guz ii- 
it follows that for t £ [0, T] — B (where B is a subset of [0, T] of 0-Lebesgue 
measure), g(.,t) £ L p { fi) and dtu(.,t ) £ L p (fl). We deduce from theorem 
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12.21 and the argumentation related to equation (12.221) that on [0, T] — B, 
u o F^ 1 (.,t) G Wp P (Q) and 


uoF (.,i)|| 


C, 


P < 


n,p,fl || 


I P 

I L°°(fl T ) 




mu(..m r L , (n) 


Integrating with respect to time we obtain that 

uoF- 1 £ LP( 0, T, wj^n)) and 


+ lli/(-,f)llk (n ))- 

(2.25) 


I 77 —111 / || || L°° /||o II i II II i 

1^°-^ Hlp(o,t,vi/? ) ’ p (o)) — 71 (l|o't u llLP(n T ) + ||9||LP(n T )J- 


1-/3, 


(2.26) 

It remains to show that under the assumptions of theorem O d t u G 
LP(Q t ). 

In order to bound ||/?tu(., i)|| LP ^ we use general Sobolev inequalities 
(chapter 5.6 of j3H]). 

• If n > 3, write p* = 2 n/(n — 2). We have for 2 < p < p*, 

( [ (d t u) p dx) § < C7 n ,n( / («9t<* dx)^ (2.27) 


thus, using Gagliardo-Nirenberg-Sobolev inequality 

1 


(< 9 t u) p dx) p < C niPi n 


-\nin(®) 


a[d t u]. 


(2.28) 


If n = 2, we write for (xi, X2, X3) G fix (0,1), x(xi, X2, X3) := dftt(xi, X2). 
Using Gagliardo-Nirenberg-Sobolev inequality in dimension three we 
obtain that for 2 < p < 6 


(d t u) p dx)p < C np & / (Vd t u) 2 dx. 


(2.29) 


Which leads us to (l2~28l) . 

If n = 1 then using Morrey’s inequality we obtain that with 7 := 1/2, 

1 


c°,T(n) 




^min(®) 


a[d t u}. 


(2.30) 


We conclude the proof of theorem 11.31 by using lemma 12.31 

We deduce theorem II .41 from Morrey’s inequality and theorem II .31 
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Holder continuity for n > 3 or non-convexity of fL In this 
paragraph we will not assume to be convex. Let N P,X (D) (l<p<oo,0< 
A < n) be the weighted Morrey space formed by the functions v : —* R 
such that ||t;|| jVP.A(n) < 00 with 

(2.31) 

To obtain the Holder continuity of u o in dimension n > 3 we use 
corollary 4.1 of m- We give the result of S. Leonardi below in a form 
adapted to our context. Consider the Dirichlet problem (12121) . We do not 
assume H to be bounded. We write W 2,P,X {D) the functions in W^ P (Q) such 
that their second order derivatives belong to N P,X (D). 

Theorem 2.3. There exist a constant C* = C*(n,p, A, <9fl) > 0 such that if 
/3 m < C* and f £ N P,X (Q) then the Dirichlet problem 12. DA) has a unique 
solution in W 2,p,x n wl' p (D). Moreover, ifO<X<n<p then X7v £ C a (Q) 
with a = 1 — n/p and 

l|Vv||c7“(n) < A ■ (M) (2.32) 

where C = C(n,p, A, (Ml). 

The proof of theorem 11.51 is an application of theorem 12.31 We just need 
to observe that from Holder inequality we have for 0 < e < 0.5 

ll/llvp> e (n) 5; C , n,p,r2,e||/||z / p(i+e)(Q)- (2.33) 

From this point the proof is similar to the proof of theorem II .31 


l w lljvp>*(n) := SU P 

xo£& 


\x — Xq 


-A 


v(x)\ 1 


2.1.2 Medium with a continuum of time scales. 

We will need theorems 1.6.2 and 1.6.3 of M- For the sake of completeness 
we will remind those theorems below in version adapted to our framework. 
Consider the following parabolic problem: 


d t v = ^2 M i:j (x)didjV + f. (2.34) 

* 4=1 


We assume M to be symmetric bounded and elliptic and v = 0 at t = 0 and 
on the boundary <9fL Write 


Trace l 1 MM] + 1 

T]m ■= sup ---o- 

xeci T (Trace[Af] + l) 


(2.35) 
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and 


Trace [M] + 1 

aM ' xeX Trace [*MM] + 1' 


Write for p > 2 

s P (n T ) := {u g LP(o,r, wl*{nj)-,dtv g l p ( 12 t );u(.,o) = 0 } 


and 


'Sp(n T ) 


' SIt 


' didjv) 2 + (<%u) 2 ) 2 dy dt. 




(2.36) 


(2.37) 

(2.38) 


Theorem 2.4. Assume 12 to 6e convex and that there exists e G (0,1) such 
that pm < l/(n+e), then for each f G L 2 (Qt ) tde Cauchy-Dirichlet problem 
\2.3f\) admits a unique solution in S^I^t) which satisfies the bound 

IMIs 2 (fi T ) < pll/llL 2 (n T )- (2.39) 


Theorem 2.5. Assume 12 to be convex and that there exists e G (0, 1) such 
that pm < l/(n + e), then there exists a number po > 2 depending on 12,n,e 
such that for each f G L p (Dt) the Cauchy-Dirichlet problem \2.32\) admits 
a unique solution in S p (Dt) which satisfies the bound 

IMIs p (n T ) - _ ^— f \\f\\LP(Q T )- (2.40) 


Remark 2.2. In fact theorem 1.6.3 of m is written with 1 — C(p)y /1 — e 
in the denominator of (non but it is easy to modify it to obtain (E35I) by 
lowering the value of po changing the value of C p . 

Let 6 > 0. Let us now apply theorem 12.41 on [0,T/<5] with 

M \= 5a o F~ 1 (y, dt) (2.41) 


and 

/ := ${g o F~ 1 )(y, 5t). (2.42) 

Observe that if condition o is satisfied then F is an homeomorphism and 
M is well defined, bounded and elliptic. Moreover t)m < oo and < oo 
since 


esssupQ T 

T 


Trace [*MM] + 1 
(Trace [M\ + l) 2 


= esssup nr 


5 2 Trace ^cru] + 1 
(<5 Trace[cr] + l)“ 


(2.43) 


It follows that the following equation admits a unique solution in 52(12 t). 


d t w(y,t) =^2M iJ (y,t)d i d j w(y,t) +k(y,t) 
id 


(2.44) 
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with k(y,t ) = 5g(y,5t). And we have 


[ S [ {{dtwf +^2(didjw) 2 ) dydt < — - C - ||/||^(a r )~ (2-45) 

Jo Jn “J {L-y/L-ey t 

Using the change of variables t —> 5t and writing 

w(y,t) ■■= v(y,5t). (2.46) 

we obtain that v satisfies the following equation on 

dtv(y, t) =Y^{v( F ~ 1 {y,t),t)) iJ d i d j v(y,t) + g(y,t). (2.47) 

Using the change of variable y = F{x) and observing that dtF = div aVF 
we obtain that v o F satisfies 

dt(v o F) = div (aV(u o F)) + g. (2.48) 

It follows from the uniqueness of the solution of that u = v o F. 

In resume we have obtained theorem n~Hi (we use lemma vn\ to control the 
constants). The proof of II.71 is similar and based on theorem 12.51 The proof 
of II .81 follows from rm and Morrey’s inequality. 

Let us now prove proDosition ll.il Write x = Trace [cr] and z = n 
(observe that 1 < z < n). It is easy to check that condition II .1 1 can be 
written 

— S 2 x 2 ( e H z — 1) + 2x5 — (n + e — 1) > 0. (2.49) 

n 

Choose 5 = n||(Trace[cr]) -1 || Loo ^^. Observing that 5x > n and 5x < ny a 
it is easy to conclude the proof of proposition 11.11 Similarly obtains the 
following lemma by straightforward computation from equation drm 

Lemma 2.4. Assume that condition [T~7] is satisfied then < C(n,e,6) 

||(Trace[o-]) _1 || LOO(OT) < C[n,e,S) (2.50) 

and 

|| Traced]|| LOO(nT) < C(n,e,S). (2.51) 

2.2 Convergence of the finite element method. 

Write 1Zh the projection operator mapping L 2 (0, T; Hq (0)) onto \fi defined 
by: for all v G Yt 

At[v,u — lZhu\ = 0. (2.52) 

Write p := u — TZhU and 9 := IZ^u — Uh- 
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Lemma 2.5. 


-\\(u-u h )(T)\\ 2 L2 ^+AT[u-u h }= / pd t (u-u h )+A T [p,u-u h \. (2.53) 
" «/ 

Proof. Subtracting (11.111 (integrated against if) and (11.311) we obtain that 


{if, d t {u - u h )) + a[if, u-u h ]= 0 

for all if G Vh{t). 

(2.54) 

Integrating by parts with respect to time we deduce that 


{if, (u - u h )(.,t)) + a[if, u - u h ] 

= / d t if{u — Uh). 

J n t 

(2.55) 

Taking if = 6 in (12.551) we deduce that 




|| {u - u h ){., t) || +A t [u -u h }= d t Q(u - u h ) + {p, (u - u h )(.,t)) 

+ A t [p,u-u h \. 

(2.56) 


Observing that 

{d t 6,u-u h ) + {p, (u-u h )(.,t)) 

we deduce the lemma. 



/•t 

(u-u h )(.,t)||* 2(n) +J 


(p,d t {u-u h )). 

(2.57) 

□ 


2.2.1 Time independent medium. 
Lemma 2.6. 


(u 


« fc )(r)||= a(n) 


+ At[u — Uh] < 2 


(l|p||L2(n T )ll^ - d t u h || L 2 ( n T ) 

+ At[p\Y 


(2.58) 


Proof. Lemma f2.6l is a straightforward consequence of lemma l231 and Cauchy- 
Schwartz and Minkowski inequalities. □ 

Lemma 2.7. We have 


IM'.rlllhn) +*W < T^fMsll| w (2.59) 

/ 'min \ u ) 

Proof. Taking if = u,h in (HSJ) and integrating with respect to time we 
obtain that 


^!K(-’ T )|Il 2 {n) +A T [u h \ = {u h ,g) L 2 (n T ). (2.60) 

Using Poincare and Minkowksi inequalities leads us to ESP - □ 
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Lemma 2.8. Assume dta = 0. We have 


\^t U h\\ L 2^ T ^+ a [uh(-,T)] < |MIl2(q t )- 


(2.61) 


Proof. The proof is similar to lemma 12.21 We need to take if = dtuy, in 

(OTT) . □ 


/ 

Jn 


Let t G [0, T] and v G Hq(Q,), we will write !Zh,w{-,t) the solution of: 

t ’Vifa(x, v — TZh^v) dx = 0 for all if G V/,.(f). (2.62) 

We will need the following lemma, 

Lemma 2.9. Assume the mapping x — » F(x,t ) to be invertible, then for 
v G Hq(Q) we have 

• For n = 1, 

(a[u - n hjt v}) 3 < C'x/i||'yoF“ 1 (., ^11^2,2 HaVFlIloo^)- (2.63) 

• For n > 2, 

(a[v - U hi tv\y <C x h\\v o F^{.,t)\\ w 2,2 

n —1 re —2 

xCnfia 4 ||(TraceH) _1 || L ^ (nT) . 


(2.64) 


Remark 2.3. Recall that p a is given by equation (II. 121) and it is easy to check 
that is bounded by an increasing function of (1 — /So-) -1 . 

Proof. Using the change of coordinates y = F(x, t) we obtain that (we write 


v := v o F x ) 

a[v] = Q[v\ 

(2.65) 

with 

Q[w\ := [ t X7w(y,t)Q(y,t)X7w(y,t)dy 

Jn 

(2.66) 

and 

Q(!/J> := det(VF) ° F “• 

(2.67) 

Using the definition of 1Zht v we obtain that 

Q[v - U h t v\ = inf Q[v - p}. 

(2.68) 

Using property (11.261) we obtain that 

Q[v~ K h ,tv] < X max (Q)Cxh 2 \\v\\ 2 w 2 , 2 {T) . 

(2.69) 

It is easy to obtain that 
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• n = 1. 

Amax(Q) < 110-V11 ^oo (n r )- (2.70) 

• n > 2. 

n — 1 n 1 

Amax(Q) < C n p a 2 \| (Trace[cr]) _1 112cx>(Q T )• (2.71) 

□ 

Lemma 2.10. Assume thatdta = 0, tt is convex, (3 a < 1 and (Trace [cr]) _1 £ 
then 

A T [p\<Ch 2 \\g\\ 2 LHnT) . (2.72) 

Remark 2.4. The constant C depends on Cx, n, 11, A m i n (a) and 

|| (Trace[< t]) _ 1 || ioo ^ r ). If n > 5 it also depends on || Trace [cx]|| LOO , f ^ and if 

n = 1 it also depends on A max (o). 

Proof. The proof is a straightforward application of lemma, Inl and theorem 
o Observe that in dimension one aS7F = (f^a x ) 1 □ 

Lemma 2.11. Assume thatdta = 0 , ft is convex, f3 a < 1 and (Trace [<t]) _1 £ 
L°°(Qt) then 

l!plU 2 (n T ) < C'^ 2 ||5 r llL 2 (o T ). (2.73) 

Remark 2.5. The constant C depends on Cx, n, 0, A m i n (a) and 

|| (Trace[a]) -1 || i00 ( n ^. If n > 5 it also depends on || Traceand if 

n = 1 it also depends on A max (a)- 

Proof. The proof follows from standard duality techniques (see for instance 
theorem 5.7.6 of fTTjh We choose v £ L 2 (0, T, 77q( 12)) to be the solution of 
the following linear problem: for all w £ L 2 (0, T, Hq (12)) 


A T [ W ,v] = (w,p) L 2 {nT y 

Choosing w = p in equation (12.7411 we deduce that 

(2.74) 

d l 2 (n T )) =A t [p,v- R h , t v\. 

Using Cauchy Schwartz inequality we deduce that 

(2.75) 

Z 5 i 2 (n T ) ^ {At\p\) 2 {A t [v - Tl h}t v ]) 2 . 

Using theorem 11 .1 1 we obtain that 

(2.76) 

^ L 2 (0,T,w£’ 2 (n)) — C\\p Ll 2 (H t )- 
Using lemma, l2~9l we obtain that 

(2.77) 

(At[v - Rh,M) 2 < Ch\\p\\ L 2(n T y 

(2.78) 
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It follows that 


\\p\\ L 2 { n T) <Ch{A T [p}) h . (2.79) 

We deduce the lemma by applying lemma 12.101 to bound At\p\. □ 

Theorem 2.6. Assume thatdta = 0 , Q is convex, (3 a < 1 and (Trace[cr]) _1 £ 
then 

|| (u - u h )(.,T)\\ L2{Q) + ||u - ^|| L 2 ( 0jT . ff i( n) ) < Cti \\g\\ L 2 (n T y (2.80) 

Remark 2.6. The constant C depends on Cx, n, f i, A m i n (a) and 
11(Trace[cr]) —1 If n > 5 it also depends on || Trace [<t]|| loo ^^ and if 

n = 1 it also depends on A max (a). 

Proof. The proof is a straightforward application of lemmas rrmrnoiiTHi 
12.61 and 12.21 □ 


2.2.2 Medium with a continuum of time scales. 

In this subsection we will assume that the finite elements are in H 2 (Q) n 
Hq (fl) and satisfy inverse inequality (11.271) . 

Lemma 2.12. 


1 , 


(u - u h )(t )||^ 2(n) + A t [u - u h } 


ln t I det VF| o F~ l 


{g + ^2 <?i,j o F~ 1 d i djU h - d t u h ). 
i,j =1 

(2.81) 


Proof. Consider equation (12331) . We have 


[ pd t {u-u h ) = f p. — dt(u-u h )+ [ pd t F{VF) l V{u-u h ). 

.hit Jsh | det VF\ of 1 J Qt 

(2.82) 

Using equation o we obtain that 


' Qt 


n » 

pd t F(yF)~ l \7(u - u h ) = —A t [p,u — u h \ - Y] / pQ hJ 8^(11 - u h ). 


(2.83) 

□ 


Lemma 2.13. 

11 d t u h 


" | det(VF)|i o F - 1 “ 2 IMIi 2 (n T ) + c \\ u h\\ L 2^ T)W ^ m y ( 2 - 84 ) 

where the constant C depends on n, A max (a), || Tracefu]||L°°(n T p p a - 
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Proof. Using the change of variable y = F(x, t) in (II .31 1) we obtain that for 
all p G Xji 

(75 detiVF^oF- 1 )l 2 (P.) ~ Si,j=l Qip^i^j^h) L 2 (Q.) 

' +(^> |det(VF)|oF- 1 )l 2 (CI) (2.85) 

U h (x, 0) = 0. 


Recall that Q is given by G23- We choose p = dtu 

a a i 

-j- =--1 det a \ 4. 

|detV.F |2 |det<r|4 


and observe that 

( 2 . 86 ) 


Thus 


II i || ^ C (fl) \ max(^)? II Trace[cr]|| L oo (f2 T ),p a )- (2.87) 

| det VF| 2 

We deduce the lemma by Minkowski inequality. □ 

Combining lemma 12.121 and lemma 12.151 we obtain the following lemma 

Lemma 2.14. 

2 IK U — Uh)(T)\\ 2 L2m + A t [u - u h ] <||p|| L 2 ( n T) (H s'Hlrjv) 

( 2 . 88 ) 

where the constant C depends on n, A max (a), || Trace[cr] || i0 om T )> Pa- 
Lemma 2.15. Assume that is convex, and condition o is satisfied then 


WpWl^Qt) < Ch 2 \\g\\ L 2^ nT y (2.89) 

Remark 2.7. The constant C depends on Cx, n, fl, 5 and e, A min (a) and 
Amax(fl) • 

Proof. The proof is similar to the proof of 1emma l2.1 11 As in 12.1 ll we choose 
v G L 2 (0, T, Hq (fi)) to be the solution of the following linear problem: for 
all w G L 2 (0, T, Hq(Q)) 

At[w,v] = (w,p) L 2 (slT) . (2.90) 

Choosing w = p in equation we deduce that 

\\p\\h(nr)) = A t[Pi v - Kh, t v]. (2.91) 

Using Cauchy Schwartz inequality we deduce that 

l|p|li 2 (n T ) < {AtIpDHMv - K h , t v\y. (2.92) 
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Using theorem 11 .61 we obtain that 


II^IIl 2 (o,t,w^ 2 (Q)) — C , |l/ > lli 2 (n T )- (2.93) 

Using lemma 127)1 we obtain that 

(At[v - 'R-hjv }) 2 < C'/i||p|| L 2( 0t ). (2.94) 

It follows that 

IHU 2 (n) <Ch(A T \p\)*. (2.95) 

We deduce the lemma by applying lemma 12.91 and theorem 11.61 to bound 
At\p\. □ 

Lemma 2.16. Assume that U is convex, and that Tracefa] £ 

C 

ll^llL 2 (o,T,w/ 2 ’ 2 (n)) — l 2 (ci t )- (2.96) 

Remark 2.8. The constant C depends on Cx, n, fl, A m i n (a) and A max (o) 
and || Trace[cj]|| l°o(O t )- 

Proof. Using the inverse inequality 071) of the finite elements we obtain 
that 

C x 

II^IIl 2 (o,t,w 2 ’ 2 (Q)) — _ /) _ II^ 7 '“ ,1 IIl 2 (o,t,tv 2 ’ 2 (o))- (2.97) 

Using the change of variables y = F(x) we obtain that 

l! v */i|li2( 0i T,vu 2 ’ 2 (n)) - c A T [uh] (2.98) 

where C depends on n, A m i n (a) and || Trace[cr]We deduce the 
lemma by using lemma 1?~71 □ 

Theorem 2.7. Assume that U is convex, and condition U . 1\ is satisfied then 

\\\(u- u h ){T) ||^ 2 (0) + A t [u - u h \ < ChM^y (2.99) 
Remark 2.9. The constant C depends on Cx, n, U, 5 and e, A min (a) and 

Amax(O'). 

Proof. The proof is a straightforward application of lemma 12.141 lemma !2.15l 
and lemma mni □ 
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2.2.3 Homogenization and discretization in time. 

We use the notation of subsection ll.3l First let us observe that the numerical 
scheme associated to (081) is stable. Indeed choosing ^ = v n+ \ one gets 

| ^n+l (tn+l) | = { v n+l(t n )> v n{tn)) ^2^ 


1 


+ 2 (l t ’ n + 1 (^ n + 1 )L 2 (n) |^n+ 1 (in)| i 2( n )) 


2 

ptn+l / \ 

/ (a[u n+ i(t)] + (v n+ i(t),g{t)) L2[n) j dt. 

J t n 


( 2 . 100 ) 


It follows by Cauchy-Schwartz inequality that 

2 I Vn+1 ( in +!) I l 2 (fi) — 2 I Vn I l 2 (H) 

- / (a[u n+ i(f)] + (u n+1 (t),ff(t)) i2(n) J dt. 

J tn 

Hence using Poincare and Minkowski inequalities one obtains that 

2 Fn+l 2 

l^n+1 (tn+l) "L I d [u n -)_l (t)] dt < |u n (f n ) 1 


( 2 . 101 ) 


+- 


a 


n,Q 


rtn+l 

/ k*)| 

J tn 


( 2 . 102 ) 


L 2 (n) 


dt. 


^min(o) 

which implies theorem II .111 and the stability of the scheme. Integrating 
on with respect to time we obtain that for if) E V, 

ptn+l / 

u h(tn+l)) = {' l l ) (t n ), Uh(t n )^ L + / ^(dt'*/’(^)> u h^)) L 2 (Q.) 

J t n 

- a[i/)(t),u h (t)\ + (i/)(t),g(t )) L2(n) ) dt. 

(2.103) 

Let us write (u 1 ) the coordinates of Uh associated to the basis (tpi o F), i.e. 

Uh(x, t ) := ^2 u\t)tpi(F(x, t)). (2.104) 


Let us define 


(2.105) 


u n (x,t) := J2^(tMF(x,t)). 

i 

Subtracting (12.1031) and Pi we obtain that for if: E Zt, 

(V’C^n+l)) ( u n +1 ^n+1) (^n+1)) = (V’(^n); (^n ^n) (tn)) 

ptn -\-1 / 

+ ((QnKf),(«h -Un+l)(*)) £ 2 (n) 

J t n 

- a[^(t), [u h - u n+ i) 


dt. 


(2.106) 
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Choosing ^ = u n+ \ — v n+ \ we deduce using Cauchy-Schwartz inequality 
that 

1 2 f tn+ 1 

- | {Un+l Ui+1) (tn+l) T / ®[(^n+l ^re+l)(t)] tit — 

■ t n 

1 2 / 

2 I (^ n — v n)(tn) |£2+ / ^(3t(u n _|_i — V n +i)(t), (llh — U n _|_i) (t)) 

i/ t n 

-a[(u n+ 1 - u n+ i)(t), (u h - u n+ i)(t)]^j dt. 

(2.107) 

Time independent medium. Observe that if the medium is time in¬ 
dependent then (12.1071) can be written 

1 2 f tn+ 1 

-|K+1 Ui+1) (tn+l) | £ 2 (Q) T / ®[(tin+l ^re+l)(t)] dt U 

tn. 


;\{u n v n){t n ) | L 2 (f 2 ) 


P^n+1 


a [(«„ + i - u n+ i ){t),(u h - u n+ 1 )(t)] dt 

(2.108) 


which leads us to 


1 2 rtn+i 

^'|( , ^n+l ^n+l) (^n+l) | “1“ J &[(^n+l ^n+l)(^)] dt ^ 


2 I t ’w)(^)| i 2( f2 ) 

rtn + 1 rt-n+l 

+ / l(t < s) a[(u n+ i - v n+1 )(t),d s u h (s)] dsdt. 

•J tn J tn 

Write At := t n +1 — t n . Using Minkowski inequality we obtain that 
a[(u n +i -v n+ i)(t),d s u h (s)\ <^^a[(u n+ i -u n+ i)(t)] 

+ -Ata[d s u h (s)}. 


(2.109) 


( 2 . 110 ) 


It follows from (12.1091) that 

|(^n+l — W+l)(tn+l)| i 2 (Q^ + 

| (Un — v n){t n ) |^ 2 (Q) + l^t| j 


t-n +1 


tn+l 


a[(u n+ i - u n+ i)(t)] dt < 
a[d s u h (s)j ds. 


( 2 . 111 ) 


Observing that 

rtn +1 rtn +1 

/ a [(u n +i - v n+ i) (t)] dt >0.5 / a[(u h - v n+ i)(t)] 

Jtn Jtn 

f>t. 


dt 


ftn +1 

/ a[{u h -u n+l )(t)\dt 

J t n 


( 2 . 112 ) 
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and 




f*tn -\-1 


a[(u h -u n+1 )(t)]dt <\At\ 2 I a[d s Uh(s)]ds (2.113) 

Jt n Jt n 

we obtain that 

2 rt-n+l 

| (^n+l ^n+l) (^n+l) "b 0.5 J ^n+l)(^)] dt — 


Cn+l 


(2.114) 


| (^n Vn){tn) \ 4” ^|At| j a[d s Ufo(s)] ds. 

J tn 

In conclusion we have obtained the following following lemma 
Lemma 2.17. Let v E &e the solution of FTM) . We have 

f T f T 

- v)(T)\\ 2 l2 (q) + / a[{u h - v)(t)\dt < 3\At \ 2 / a[9 s u ft (s)] ds. 

(2.115) 

Combining lemma, l2~8l with lemma, 12.171 we obtain the following theorem: 
Theorem 2.8. Let v E Zt be the solution of \1.3H\) . We have 

\\(uh ~ v ){T)\\ 2 L 2 {n) + a[{u h -v){t)]dt < 3|At| 2 

+ l|5(->°)||L 2 (n))- 

(2.116) 

Time dependent medium. Observe that 

dti^n+l ^Vi+l) — dtF{S7F^ \7(u n+1 V n -\-i). 

It follows after writing dtF = divaVF, integration by parts and using the 
change of variables y = F(x,t) in (12.1071) that 


2 ftn+l 

| (^n+l ^n+l)(tn+l) T / ®[(^n+l ^n+l)(t)]dt ft 

J t n 


o IK ^n)(tn)| L 2 (n ) 


£n+l 


- 2 


a[(«n+l - Wn+l)(t), (w/i - U n+ i)(t)] dt 


E 

ij ' 




{uh Hn+l)Qi,jdidj{hn+\ ^n+1 )dtdy. 


(2.117) 
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Hence using Minkowski inequality we obtain that 

2 ftn+l 

| (^n+l ^ri+1) (j'n+l ) | "h J &[(^n+l dt ^ 

2 pt n +i 

\(u n -v n )(t n )\ L2 ^ + 4 J a[(u h -u n+ i)(t)]dt 

rtn+i r 

2 / / (“U/i ^n+i )Qi,j&iOj(ii n +\ u n _(_i) dtdy. 

^ j ^ tn J ^ 

Using Minkowski inequality we obtain that 

| /* ^n+l /* 

^ ^ / / (^/i ^'n+l)Qi, < 7^i^7(^n+1 ^n+l) dtdy ^ 

• • tn J fl 


CUrc : 


/•*n+l r 

J tn «/ O 


u h - u n+ 1 | 2 dtdy 


+ 


^max (Q) 

Ca 


rtn+i _ r 

/ XU l<%< 9 j(^n+l - •O n+ i)| 2 dt(i?/. 

J tn i j J ^ 


Using the inverse inequality 071) and the change of variable y = 
obtain that 

C tn , _ _ 

< 


rtn+i _ r 

/ / \9idj{u n+ i - v n+ i)\ 2 dtdy 

^ t n t j ^ ^ 


Cx 


h 2 \ m i+Q) 


rtn+l 

/ a[(u n+ i - v n+ i)(t)]dt. 

J t n 

In resume, choosing Ca = We ^ iave obtained that 

2 rtn+i 

| (^n+1 ^n+l)(^n+l)|^2^j 0-5 / ^[(^71+1 ^n+l)(^)]^ — 

■' t n 

2 /*^n+l 

I (^n ^n) (^ 71 ) I"I - ® J 'U j n+l')(j')\ dt 


8CxA max ((2) 2 
h 2 \ min (Q) H 


rtn+i r 
J tn v O 


% - -u n+ i| 2 dtdy. 


And a computation similar to the one leading to (12.11011 gives us 

2 i rtn+i 

| (tl n +\ V n +i)(t n +i) | ~b ^ / o[(^h ^n+i )(^)] dt — 

■' tn 


2 

I (^71 ^71) (^71) | ^2(^2) H” 9 J ^71+1 )(^)] 




8CjYA max (Q) 2 

4 - - 77 

^ 2 Amin(< 3 ) 


/■*n+l r 

/ / 1 % 

w tn </ 


- u n+ i| 2 dtd?/. 


(2.118) 


(2.119) 


F(x) we 


( 2 . 120 ) 


( 2 . 121 ) 


( 2 . 122 ) 
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Moreover using the change of variables F{x) = y and the inverse inequality 
(fOKl) we obtain that 


f tn+1 U ^ / CxXmtxiQ) [ tn+1 fn . !2 jj. ; 

/ a[(u /l - u n+ i)(t)J dt < -72- / / |ith-Mn+i| dtdy. 

J t n J n 





(2.123) 

Let us also observe that 




[ tn +i r 


ftn+1 f 


/ / \u h - u n+ i\ dtdy < 

|At | 2 

1 / | d t u h \ 2 dtdy. 

(2.124) 

J tn ** ^ 

J 

tn ^ 


It follows that 




| (^n+l — v n+l){tn+l)\ L 2(ty + ^ 

r tn-\- 1 



Jtn 

rtn + 1 

a[(u h - u n+ i)(t)] dt < 

f 

(2.125) 

|K v n)(t n )\ L 2 ^ + Cb 

J t n 

/ \d t u h \ 2 dtdy. 

Jn 



with 

Cb = C' n C Y A max (Q)(l + 1 ). (2.126) 

^minAV j 


We deduce that 

\\(u h -v){T)\\ 2 L2{n) + ^ J a[(u h - v)(t)]dt < J J \d t u h \ 2 dt dy. 

(2.127) 

Using lemma, ET31 to control Cb and combining (12.1271) with theorem 11.61 we 
obtain the following theorem. 

Theorem 2.9. Assume that U is convex, and condition U.H is satisfied. Let 
v E Zt be the solution of arm . we have 

U u h ~ v)(T)\\ 2 l2{q) + a[(u h -v)(t)]dt <C^-\\g\\ 2 L2 ^ T y (2.128) 

where C depends on il, n, 6, e, A min (a) and A max (a). 
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3 Numerical Experiments 

The purpose of this section is to give several illustrations of the implemen¬ 
tation of this method. The domain is the unit square in dimension two. 
Equation ten is solved on a fine tessellation characterized by 16129 inte¬ 
rior nodes (degree of freedoms). 

Three different coarse tessellations are considered, one with 9 degrees of 
freedoms (noted dof in the tables), one with 49 and the last one with 225. 

The parabolic operator associated to equation JEU) has been homoge¬ 
nized onto these coarse meshes using the method the method presented in 
this paper. We have chosen splines to span the space Xh introduced in 
subsection PI 

3.1 Time independent examples. 

Example 1. Time independent site percolation. 

In this example we consider the site percolating medium associated to 
figure Cl The fine mesh is characterized by 16641 nodes. m has been 
homogenized to three different coarse meshes with 9, 49 and 225 interior 
nodes using the method described here and splines for the space X) } . 0) 
has been solved with the fine mesh operator and the coarse mesh operators 
with g = 1 and g = sin(2.4x — 1.8y + 2irt). The fine mesh and coarse mesh 
errors are given in tables □Email Figure □shows u computed on 16641 
interior nodes and Uh computed on 9 interior nodes in the case g = 1 at 
time 1. 



(a) u. (b) u h . 

Figure 7: u computed on 16641 interior nodes and Uh computed on 9 interior 
nodes. 
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Tabic 1: Coarse A/ 


esh Error. Time Independent Site Percolation with g — 1. 


wsm 

L 1 

L°° 

L 2 

H 1 

9 

0.0142 

0.0389 

0.0168 

0.0366 

49 

0.0077 

0.0450 

0.0101 

0.0482 

225 

0.0035 

0.0228 

0.0060 

0.0293 


Table 2: Fine Mesh Error. Time Independent Site Percolation with g — 1. 


Mtsm 

L 1 

L °o 

L 2 

H 1 

9 

0.0196 

0.0843 

0.0251 

0.1193 

49 

0.0136 

0.0698 

0.0184 

0.1028 

225 

0.0040 

0.0243 

0.0070 

0.0485 


Table 3: Coarse Mesh Error. Time Independent Percolation Case with g 
sin(2.4x — 1.8y + 2nt). 



L 1 

L°° 

L 2 

H 1 

9 

0.0236 

0.0569 

0.0262 

0.0477 

49 

0.0181 

0.0571 

0.0215 

0.0558 

225 

0.0119 

0.0774 

0.0167 

0.0939 


Table 4: Fine Mesh Error. Time Independent Percolation with g = sin(2.4x 
1.8y + 2nt). 


MUM 

L 1 

L°° 

L 2 

H 1 

9 

0.0424 

0.1099 

0.0512 

0.1712 

49 

0.0277 

0.0985 

0.0348 

0.1451 

225 

0.0174 

0.0886 

0.0242 

0.1192 
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Table 5: Coarse Mesh Error, high conductivity channel. 


dof 

L 1 

L°° 

L 2 

H 1 

9 

0.0159 

0.0496 

0.0207 

0.0477 

49 

0.0067 

0.0389 

0.0102 

0.0345 

225 

0.0035 

0.0228 

0.0060 

0.0293 


Table 6: Fine Mesh Error, high conductivity channel. 


dof 

L 1 

L°° 

L 2 

H 1 

9 

0.0178 

0.0564 

0.0257 

0.0947 

49 

0.0079 

0.0388 

0.0129 

0.0660 

225 

0.0040 

0.0243 

0.0070 

0.0485 


Example 2. Time independent high conductivity channel. 

In this example a is random and characterized by a fine and long ranged 
high conductivity channel. We choose a(x) = 100, if x is in the channel, and 
a(x) = 0(1) and random, if x is not in the channel. The media is illustrated 
in figure 0 



Figure 8: High Conductivity Channel superposed on a random medium. 


Tables El and El give the coarse and fine meshes errors. 
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3.2 Time dependent examples. 

In the following examples we consider media characterized by a continuum 
of time scales. In the following examples the ODE obtained on the coarse 
mesh from the homogenization of O) have also been homogenized in time 
according to the method described in subsection 11,31 

Example 3. Time Dependent Multiscale trigonometric. 

In this example a is given by equation cm Although the number fine 
time steps to solve o> is 2663, only 134 coarse time steps have been used 
to solve the homogenized equation. Hence if one also takes into account 
homogenization in space, the compression factor is of the order of 35000 for 
the coarse mesh with 9 interior nodes. 

Figure El shows the curves of t —» a(xo,t) and t —► F(xo,t ) for a given 
xq G P. 




(a) t -*■ a(x 0 ,t). 


(b) Top view of t 

F(x 0 ,t). 


Figure 9: Multiscale time dependent trigonometric medium. 


The coarse and fine mesh relative L 1 , L 2 , L°°, and H 1 errors with respect 
to time have been plotted in figures II011111II21 and II31 The initial increase 
of the relative error has its origin in the initial value u = 0 at time 0. 

The coarse and fine meshes errors are given in tables0and|H]for g = 1 at 
t = 0.1, those errors are given in tables IHland 1 101 for g = sin(2.4x — 1.8y + 27rt) 
at t = 0.1 
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Table 7: Coarse Mesh Error. Multiscale trigonometric time dependent Medium. 
<7 = 1- _ 


dof 

L 1 

L°° 

L 2 

H 1 

9 

0.0018 

0.0045 

0.0019 

0.0039 

49 

0.0012 

0.0054 

0.0015 

0.0060 


Table 8: Fine Mesh Error. Multiscale trigonometric time dependent Medium. 
9 = !• _ 


dof 

L 1 

L°° 

L 2 

H l 

9 

0.0031 

0.0096 

0.0034 

0.0242 

49 

0.0014 

0.0059 

0.0016 

0.0166 


Table 9: Coarse mesh error. Multiscale trigonometric time dependent Medium. 


- 27 rt] 

. 

dof 

L 1 

L°° 

L 2 

H 1 

9 

0.0043 

0.0087 

0.0044 

0.0085 

49 

0.0033 

0.0079 

0.0035 

0.0084 


Table 10: Fine mesh error. 
g = sin(2.4x — 1.; 


Multiscale trigonometric time dependent medium. 


- 27 rt] 

. 

dof 

L 1 

L°° 

L 2 

H 1 

9 

0.0082 

0.0199 

0.0087 

0.0379 

49 

0.0038 

0.0104 

0.0040 

0.0244 
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(a) Coarse mesh L 1 error. (b) Fine mesh L 1 error. 

Figure 10: L 1 error. Multiscale Trigonometric time dependent Medium. 




(a) Coarse mesh L 2 error. (b) Fine Mesh L 2 error. 

Figure 11: L 2 error. Multiscale Trigonometric time dependent Medium. 
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coarse mesh L 




(a) Coarse mesh Ler- (b) Fine mesh error, 

ror. 

Figure 12: L°° error. Multiscale Trigonometric time dependent Medium. 


mesh H t error fine mesh H ( error 




(a) Coarse Mesh // ' er- (b) Fine Mesh // 1 error, 

ror. 

Figure 13: H 1 error. Multiscale Trigonometric time dependent Medium. 
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Example 4. Time Dependent Random Fourier Modes. 

In this example a(x, y, t ) = e h ( x,y,t ^ where h is given by the following 
equation 

h(x,y,t ) = E (afc sin(27rfc.x / ) + bk cos(2irk.x')) 

\k\<R 

where R = 4, x' = x + \/2f, y' = y — \/2i, Ofc and are independent 
identically distributed random variables on [—0.2,0.2], In this example, one 
can compute that * max f a ) = 95.7. 



Figure 14: Random Fourier Mode 


Figure mi is a plot of a at time 0. Although the number fine time steps 
to solve m is 1332 only 67 coarse time steps have been used to solve the 
homogenized equation. Hence if one also takes into account homogenization 
in space, the compression factor is of the order of 35000 for the coarse mesh 
with 9 interior nodes. 

The coarse and fine mesh relative L , L 2 , L°°, and H 1 errors with respect 
to time (up to time t = 0.1) have been plotted in figures H5l 17771H7I and H8l 
Those errors are also given up to t = 1 in figures 1191 12U117111771 

The coarse and fine meshes errors are given in tables and El for g = 1 

at t = 0.1. 
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(a) coarse mesh L 1 error. (b) fine mesh L l error. 

Figure 15: L 1 error. Random Fourier modes up to t — .1. 


Table 11: Coarse mesh error at t — 0.1. Time dependent random Fourier modes. 


dof 

L 1 

L°° 

L 2 

H 1 

9 

0.0025 

0.0028 

0.0064 

0.0052 

49 

0.0032 

0.0098 

0.0034 

0.0100 


Table 12: Fine mesh error at t = 0.1. Time dependent random Fourier modes. 


dof 

L l 

L°° 

L 2 

H 1 

9 

0.0063 

0.0344 

0.0079 

0.0481 

49 

0.0042 

0.0207 

0.0049 

0.0337 
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coarse mesh L error fine mesh L error 




(a) coarse mesh L 2 error. (b) Fine Mesh L 2 error. 

Figure 16: L 2 error. Random Fourier modes up to t — .1.. 




(a) coarse mesh error. (b) fine mesh / JOC error. 

Figure 17: L°° error. Random Fourier modes up to t — .1.. 
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fine mesh H 






(a) coarse mesh H 1 error. (b) tine mesh H 1 error. 

Figure 18: H 1 error. Random Fourier modes up to t = .1.. 




(a) coarse mesh L 1 error. (b) tine mesh L 1 error. 

Figure 19: L 1 error. Random Fourier modes up to t — 1. 
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coarse mesh L error fine mesh L error 




(a) coarse mesh L 2 error. (b) Fine Mesh L 2 error. 

Figure 20: L 2 error. Random Fourier modes up to t — 1. 




(a) coarse mesh / JOC error. (b) fine mesh error. 

Figure 21: L°° error. Random Fourier modes up to t — 1. 
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coarse mesh H f error fine mesh H ( error 




(a) coarse Mesh H 1 error. 


(b) fine Mesh H 1 error. 


Figure 22: H 1 error. Random Fourier modes up to t — 1. 
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Tabic 13: Coarse Mesh Error for the time dependent random fractal medium with 
spline element 


dof 

L 1 

L°° 

L 2 

H 1 

9 

0.0046 

0.0074 

0.0052 

0.0065 

49 

0.0036 

0.0046 

0.0036 

0.0059 


Table 14: Fine Mesh Error for the time dependent random fractal medium with 
spline element 


dof 

L 1 

L°° 

L 2 

H 1 

9 

0.0039 

0.0082 

0.0043 

0.0222 

49 

0.0033 

0.0054 

0.0034 

0.0168 


Example 5. Time Dependent Random Fractal 

In this example, a is given by a product of discontinuous functions os¬ 
cillating randomly at different space and time scales. Namely a(x, t ) := 
a±(x, t)a 2 (x, t) ■ ■ ■ a n (x, t ) with n = 6 and aj(x, t) = c l pq (t ) for x £ [Jy, x 
[~, 2il). The coefficients c l pq (t ) are chosen at random with uniform law in 
[^, 7 ] with 7 = 0.7 and independently in subdivision in space and in time, 
thus they are assumed to be constant in each time interval 0.1 x [■—, ^jr). 
In this example we have = 160.3295. Although the number of fine 

time steps to solve o> is 3482, only 175 coarse time steps have been used 
to solve the homogenized equation which corresponds to a reduction of the 
complexity of the scheme by a factor of 35000 in the case of the coarse tes¬ 
sellation with 9 interior nodes, o and the map (F\, F?) are drawn in figure 
E 2 L 1 , Li 2 , Lqo and H\ errors are given in figure 1^71 to EH Coarse and fine 
mesh errors are given in table ED and El at time t = 0.1. We have chosen 
g = 1 in this numerical experiment, one obtains similar results by choosing 
g = sin(2.4x’ — 1.8y + 2nt). 
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(c) a at t = 0.1. (d) (F-\ , F?) at t = 0.1. 

Figure 23: a and (F\, F 2 ) at time t = 0, t — 0.1 for the time dependent random 
fractal medium. 
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(a) coarse mesh L l error. 


(b) fine mesh L l error. 


Figure 24: L 1 error for the time dependent random fractal medium at t = 0.1. 




(a) coarse mesh L 2 error. 


(b) Fine Mesh L 2 error. 


Figure 25: L 2 error for the time dependent random fractal medium at t — 0.1.. 
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line mesh L 




(a) coarse mesh / JOC error. (b) fine mesh / JOC error. 

Figure 26: L°° error for the time dependent random fractal medium at t = 0.1. 




(a) coarse Mesh HI error. (b) fine Mesh H 1 error. 

Figure 27: H 1 error for the time dependent random fractal medium at t = 0.1. 
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